UCSD/PTH 95-11 



Model-Independent Determinations of 
B — ► DID* D*W Form Factors 



C. Glenn Boyd*, Benjamin Grinsteim and Richard F. Lebed* 

Department of Physics 
University Of California, San Diego 
La Jolla, California 92093-0319 

We present nonperturbative, model-independent parametrizations of the individual 
QCD form factors relevant to B — > D*lv and B — > Dlu decays. These results follow 
from dispersion relations and analyticity, without recourse to heavy quark symmetry. To 
describe a form factor with two percent accuracy, three parameters are necessary, one of 
which is its normalization at zero recoil, .F(l). We combine the individual form factors 
using heavy quark symmetry to extract values for the product from B — > D*W 

data with negligible extrapolation uncertainty. 



August 1995 

* gboyd@ucsd.edu 

' bgrinstein@ucsd.edu 

* rlebed@ucsd.edu 



1. Introduction 

A determination of the Cabibbo-Kobayashi-Maskawa element \V c b\ from the 
decays B — > D*lv and B — > Dlv requires knowledge of the transition amplitudes 
(D*(p', e)|(V^ - A fJL )\B(p)) and {D{p')\V lJb \B{p)) , respectively In the limit of infinitely 
heavy b and c quark masses these amplitudes are predicted |l|-|3| at one kinematic point, 
namely, when the recoiling D* or D is at rest in the rest frame of the decaying B. In terms 
of q = p — p', this zero recoil point occurs at q 2 = g^ ax = (M — m) 2 , where M = tub 
is the decaying meson mass and m = vein or ran* is the mass of the final state meson. 
Unfortunately, the differential decay width dT/dq 2 vanishes at <7m ax , so an extraction of 
\V c b\ requires the extrapolation of the matrix element from q 2 values less than q 2 ^.^. 

This method has been used by several experiments |§-|6[]. ARGUS tested the impor- 
tance of the extrapolation on the determination of |V C &| by using various parametrizations. 
The observed variations of \V c b\ were larger than the rest of the errors combined. In prin- 
ciple the error inherent in the extrapolation can be made arbitrarily small by collecting an 
arbitrarily large amount of data, arbitrarily close to q 2 ^; such an endeavor is impractical, 
if not unattainable. Therefore, a precision measurement of |V C b| from B — > Dlv or D*lv 
requires a model-independent understanding of the extrapolation. 

In a previous letter we presented such a model-independent extrapolation. To this 
effect we used analyticity, crossing symmetry, and QCD dispersion relations to find a two- 
parameter fit to the S-meson 6-number elastic form factor F(q 2 ). Heavy quark symmetries 
were then invoked to relate F(q 2 ) to the amplitude for B — > D*W; given the validity of 
heavy quark symmetries, we showed that over the relevant range of q 2 the accuracy of the 
two-parameter fit was better than 1%. 

With a two-parameter extrapolation at hand, experimentalists can accurately deter- 
mine \V c b\ by making a simultaneous fit of the data to \V C },\ and the two parameters in our 
extrapolation. At the moment this program suffers from two main theoretical uncertain- 
ties: 

1 ) Incalculable nonperturbative corrections to the amplitudes for B — > Dlv and D*lv at 

"?max appear at orders l/m c and 1/m 2 , respectively ||. The size of these is controver- 
sial. 

2 ) The extrapolation of Ref. relies on heavy quark spin and flavor symmetries, with 

a priori corrections of order l/m c . 

Of these issues only the latter is addressed in this paper. Instead of assuming particular 
numerical values for the normalization of form factors at zero recoil and making a fit to 
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two parameters plus \V c b\, we evade the first issue by presenting our results as three- 
parameter fits in units of .F(l), the amplitude at q^ ax - We then improve on the method 
of Ref. 0] by dropping the unnecessary use of heavy quark symmetries. To this end we 
derive, in Sec. 2, bounds on the form factors describing the amplitudes for B — > Dlv 
and D*lv. As in Ref. [[7], our arguments are based on QCD dispersion relations, crossing 
symmetry, and analyticity. The bounds take the form of integrals over the unphysical 
region q 2 > (M + m) 2 , which are then related to the individual form factors in the physical 
region < q 2 < g^ ax . In Sec. 3 we derive our parametrizations by constructing quantities 
from each form factor that can be legitimately expressed as Taylor series with bounded 
coefficients. These parametrizations constitute our main results. Some technical issues are 
addressed in Sec. 4, where we demonstrate that the error incurred by ignoring cuts in the 
form factors is negligible. Section 5 enumerates the corrections to the parametrization and 
estimates their effects on the bounds. We discuss the sensitivity of our method to such 
corrections and demonstrate that their effects are minimal, thus establishing the robustness 
of the technique. In Sec. 6 we use heavy quark symmetry to relate the separate form factors 
appearing in the measured rate, but point out that measurements of individual form factors 
in the near future will obviate the need for this use of heavy quark symmetries. We present 
results from fits of current data to this parametrization, including values for | V^l^^), in 
Sec. 7. Our conclusions appear in Sec. 8. 



2. Dispersion Relations 

The QCD matrix elements governing the semileptonic decays B — > D*lv and B — ► Dlv 
may be expressed in terms of the form factors 

(D*(p\e)\V»\B(p)) = ige^e* a p' pPj 

(D* (f/, e) | \B(p)) = f e*» + (e* • p) [a+(p + p')" + a_ (p - p') M ] (2.1) 
(D(p')\ V» \B(p)) =f + (p + p'r + f-(p-p'T 

where = c'j^b, and A 11 = 07^756. In terms of these form factors, the differential decay 
widths for B — > Dlv and B — > D*lv are respectively 

dT _ \V ch \ 2 G 2 F {k 2 q 2 )i 2 

dq 2 24tt 3 M 3 IJ+I { ' ' 
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and 

dT _ |K b | 2 G^vW 



[2q 2 \f Q \ 2 + \F 1 \ 2 + 2q 4 k 2 \g\ 2 ], (2.3) 



dq 2 96tt 3 M 3 
where 

1 r _ _ 1 _ _ 1 

(2.4) 



m 



2q 2 k 2 a+-^(q 2 -M 2 + m 2 )f 



determines the partial width to longitudinally polarized D*'s, and fo and g respectively 
determine the axial and vector contributions from transversely polarized D*'s (longitudinal 
polarizations do not contribute to the vector matrix element in the B rest frame, as is 
readily seen from Eq. ( |2.1|) ). k 2 is related to the three-momentum squared p 2 D for D or 
D* in the B rest frame, and is given by 



M 2 1 

k " = -jTPd = ^1<1 2 -( M + ™?\ [Q 2 -i M - m)% (2.5) 



with M and m the B and D or D* meson masses, respectively. 

In our derivation of constraints from dispersion relations, we follow the well-known 
methods developed by authors listed in Ref. ||. We begin by considering the two-point 
function 

Uf(q) = {q^q v - q 2 g^ u )Ii T j(q 2 ) + g^U^q 2 ) = i J d 4 x e iqx (0\TJ"(x) J t!/ (0)|0), (2.6) 

where J = V or A. In QCD we can render both sides of this relation finite by making one 
subtraction. We thus obtain the once-subtracted dispersion relations 



T}L( ^ _ d\^_ = i r A imnf (i) 

dq 2 Ti J (t- q 2 



xTi<l 2 )=-^r=zl dt^-^L. (2.7) 



The functions Xj ,L (q 2 ) may be computed reliably in perturbative QCD for values of q 2 
far from the kinematic region where the current J can create resonances: specifically, 
(m& + ?b c )Aqcd <^ (Trib + m c ) 2 — q 2 . For resonances containing a heavy quark, it is 
sufficient to take q 2 = 0. 

The absorptive part ImIlj U (q 2 ) is obtained by inserting on-shell states between the 
two currents on the right-hand side of Eq. ( |2.6| ). For \i = u, this is a sum of positive-definite 
terms, so one can obtain strict inequalities by concentrating on the term with intermediate 



states of B-D or B-D* pairs. The contribution of B-D* pairs to the right-hand side of 
( |2.7|) enters (no sum on u) as 

// 7 2 
(iO % 2 - (M + m) 2 ) 

167t 2 a/q 2 ^ g N 

e 

with an analogous form (no sum over polarizations) for B-D pairs. Here rif is the num- 
ber of light valence quark flavors for the B and D or D* that give physically equivalent 
contributions; in practice, we take rif = 2. The momentum g here and subsequently is 
not to be confused with q in Eq. (|2.7|) , which will subsequently be set to zero. The matrix 
elements in Eq. ( |2.8| ) are related by crossing symmetry to those in Eq. (|2.1| ) . That is, they 
are described by the same form factors, but defined in different regions of the complex q 2 
plane, k 2 is still defined by Eq. ( |2.5| ) but may now be interpreted as the three-momentum 
squared of either the B or D, D* in the center of mass frame. For massless leptons it 
turns out that the partial widths appearing in Eq. (|2.3| ) present the same combinations of 
form factors as the space-space components of Eq. (|2.8|). It therefore suffices to use the 
dispersion relation 

Xj = -I dt f3 J{ \ (2.9) 



7T 

,T(r\\ 1 d „,L 







where Xj = Xj(0) ~ \~SpXj{ty- This definition of x,j corresponds to the combination of 
Ilj and Ilj that gives irj at q 2 = 0. At one loop, 



Xv( u ) = X A (- U ) 
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327r 2 m 2 (l 

x [(1 - m 2 )(3 + 4u- 21u 2 + 40m 3 - 21m 4 + 4m 5 + 3m 6 ) ( 2 ' 10 ) 
+ 12M 3 (2-3M + 2M 2 )hiM 2 ], 

where u = ^ is the ratio of quark masses. For u = 0.33, Xy = 9-6 • 10 _3 /m 2 and 
X^ = 5.7-10- 3 /m 2 . 

When substituted into Eq. ( |2.9|) , (|2.8| ) and ( |2.10|) lead to bounds on integrals of the 
analytically continued form factors. For example, for the axial current J = A, Eq. ( |2.8| ) 
becomes 



' - 12nJq 2 



IAI'+2>I* 
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6{q 2 - (M + m) 2 ). (2.11) 



The bound in this case, which may be taken to constrain |/q| and separately, reads 



n f 



^ X A J ( m + 



(r 



\M 2 + ^ 



< 1. 



(2.12) 



We now define a new variable z by 



1 + z 



1 - z 



(M + m) 



AMm 



(2.13) 



Taking the principal branch of the square root in this expression, the change of variables 
q 2 — > z maps the two sides of the cut q 2 > (M + m) 2 to the unit circle \z\ = 1, with 
the rest of the q 2 plane mapped to the interior of the unit circle. In particular, the real 
values — oo < q 2 < (M — m) 2 and (M — m) 2 < q 2 < (M + m) 2 are mapped to the real 
axis, 1 > z > and > z > — 1 respectively. Written in terms of z, the inequalities from 
Eqs. (ETgD- flnSD now read 



1 

2ni 



dz 



\Hz)F t {z)\ 2 <1. 



(2.14) 



c 



The contour C is the unit circle. The weighing functions are 

</>i = M 2 - s 2 2+p ^J[r(l + z)] 2 ^ (1 - z) s --z [(1 - z)(l + r) + 2^(1 + ^)]" s " p , (2.15) 



where r = m/M is the ratio of meson masses, and k, p and s depend on the form factors 
Ft as listed in Table 1. 
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1/k 
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s 
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12nM 2 x A 


1 
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Ft 


24nM 2 XA 


1 


4 


2 
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12nM 2 Xv 


3 


1 


3 


U 


6nM 2 Xv 


3 


2 



Table 1. Factors entering Eq. ( |2. 1 for the form factors Fi. 

The results (|2.14|) -( p,15|) apply equally well to analogous heavy-to-light form factors 



such as in B — > K*^ and B — > nlv; for the latter process, they agree with Ref. |I(J upon 
substitution of for mo- 



3. Parametrization of Form Factors 



Our parametrizations of the form factors rely on a Taylor expansion about z = 0. To 
connect this expansion to bounds at \z\ = 1, we need a function which is analytic inside 
the unit disk. The form factors Fi have cuts and poles along the segment q 2 > (M — m) 2 
of the real axis in the complex q 2 plane, and therefore only on the segment (—1,0) of the 
real axis in z or on the unit circle \z\ = 1. 

We have used the freedom to redefine <pi by a phase to ensure that it has no poles, 
branch cuts, or zeros in the interior of the unit circle \z\ < 1, but the form factors Fi(q 2 ) 
have poles due to the existence of stable spin-one states with unit bottom and charm 
number (spin-zero states only contribute to /_ and a_, which, for massless leptons, give 
vanishing contribution to the differential rate). The masses of these B* mesons can be 
reliably computed HTTHT3| with potential models. The vector states are predicted to have 
masses corresponding (for z defined with m = ran*) to z\ = —0.284, zi = —0.472, Z3 = 
—0.531, and Z4 = —0.907, while the axial vector masses correspond to z$ = —0.395, 
zq = —0.399, Z7 = —0.609, and z% = —0.619. One may form functions P(z) that are 
products of terms of the form (z — Zi)/(1 — Ziz), known to mathematicians as Blaschke 
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- LJ : (1 - ZjZ) 



(3.1) 



p 3 =n {z ~ Zj) 



(1 — ZjZ) 

Such P^s are analytic on the unit disk for \zj\ < 1 and serve to eliminate poles of Fi 
at each z = Zj when formed into the products Pi(z)Fi(z). Most importantly, each Pi is 
unimodular on the unit circle, and therefore we may replace Fi with PiFi in our bound 
Eq. (|2.14j ) without changing the result. Since now both PiFi and <pi are analytic on the 
unit disc, Taylor expanding ^P^Fi about z = gives 

. 00 

fl W = flf5SwS^" (3 ' 2> 

Substituting this expression into Eq. ( |2.14j ) gives the central result 

00 

$^|a n | 2 <l. (3.3) 

n=0 



The coefficients a n are different for each form factor, and must be determined by exper- 
iment. However, since both B-D* and B-D states contribute to the same vector-vector 
dispersion relation, the sum of the squared-coefficient sums for /+ and g is bounded by 
one: 

oo 

£(k /+ f + kf)| 2 )<i. (3.4) 

n=0 

This relation holds if z in Eq. (|3.2| ) is defined using m = m-^* for g and m = for / + . 
An analogous result constrains the coefficients a n of the form factors fo and Fi, 

oo 

£(k /0 ¥ + l4^| 2 )<l. (3.5) 

n=0 

For the remainder of this paper, we content ourselves with the weaker constraint Eq. ( |3.3| ). 

The utility of this parametrization arises from the observation that the physical range 
<7max > q 2 > for B — > D*lv (DW) semileptonic decays corresponds to < z < z max = 
0.056(0.065). We define an approximation F^ to the form factor Fi by truncating after 
the iVth term: 

1 N 

F l N (z)= r £a n ^. (3.6) 

Pi{z)(j)i{z) ^ 

Then the maximum error incurred by truncating after N terms is just 

^ oo 

m a x\F i (z)-F l N (z)\ = -— £ 



\Pii*)Mz)\ nr _., + i 



'a n \ z' 



< 



< 



1 z N+1 



J2 z 2n ( 3 - 7 ) 

n=N+l 



|-Pi(^max)^ , i(^max) | \/l — Z~ 



-max 



where we have used the Schwarz inequality, Eq. ( [373]) , and the fact that z N+1 /\Pi(z)4>i(z)\ 
increases monotonically over the physical range. For iV as small as 2, this truncation error 
is quite small; see Table 2. 

To calculate a relative error we need to estimate the form factor itself. This can be 
done at z = using heavy quark symmetries. The resulting bound on the relative error, 
\Fi(z) - Fl*(z) \/Fi(0), is shown in Table 2. 

The larger relative error associated with F\ arises from a collusion of factors. Com- 
pared to fo, these consist of a smaller value of n and a greater value of s +p, both of which 
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decrease <pi, as well as a smaller normalization -Fi(O). The accuracy of the parametrization 
of F\ is improved to 0.34% by truncating after one more parameter (i.e., taking N = 3 
above) . 



i 


Fi 


iF^^-Ff^lxlO 2 


\F t (z)-F z N (z)\/F(0) 





fo 


1.2 


1.0% 


1 


Fx 


4.6 


6.1% 


2 


9 


0.8 


0.5% 


3 


f+ 


1.4 


1.3% 



Table 2. Bounds on truncation errors, \Fi(z) — F i N (z)\, for N = 2, for various 
form factors from Eq. ( |3. 7f ). To estimate a corresponding relative error, we use 
the value of the form factor at threshold, Fi(0), as predicted by heavy quark 
symmetries. 



4. Branch Cuts 

In the previous section we ignored branch cuts in the form factors with branch points 
inside \z\ = 1. These cuts originate from non- resonant contributions with invariant masses 
below M + m. For example, branch points are expected at q 2 = (m^* + nm^) 2 , with n 
a positive integer, and at q 2 = (m Vbc + m p ) 2 , where rji, c is the pseudoscalar partner of the 
vector B*. We now show that their neglect is quite justified. 

We content ourselves with estimating the effect of any single cut modeled in a rea- 
sonable way, since multiple cuts can be handled analogously, and cuts modeled differently 
give comparable results. ^ 

Any form factor g(q 2 ) in Eq. ( |2.1| ) satisfies a simple dispersion relation 

i r it ijnm 

n Jo t-q z 



9(q 2 ) = ^l dt'-^-. (4.1) 



A reasonable model for a cut can be obtained, say, by taking an additive contribution to 
g satisfying 

Tmg(t) = C ^t-M 2 6(t - M 2 ) - yft - M 2 6(t - M a 2 )) , (4.2) 



' The statement in Ref. ]7j that the effect of such cuts may be incorporated by mapping them 
onto the unit circle and expanding in a new basis is erroneous; the matching of coefficients in the 
new basis to a Taylor expansion about z = involves an infinite number of equally important 
terms. 
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where q 2 = M 2 is the location of the branch point of interest, and is an arbitrary scale 
with M 2 > M 2 . The subtraction is performed to ensure that Im g(t) vanishes as t — ► oo. 
The precise form of the subtraction is immaterial, so we choose one that simplifies our 
calculations. Moreover, since branch points on \z\ = 1 are irrelevant, we need consider 
only the case M& < M + m. In all cases M a > mg» +m 7T , the location of the lowest branch 
point, which has z = —0.32 for m = mo*- 

The coefficient C arises as a coupling in diagrams connecting the (V — A) current to 
an external B-D or B-D* pair through non-resonant on-shell intermediate states. The 
intermediate states couple to the current with a strength /, and to B-D or B-D* with 
strength g. Furthermore, two-particle phase space provides a factor of 1/87T. Phenomeno- 
logically C ~ fg/8iv is expected to be quite small. We consider the most extreme case, 
namely, C = M s_3 c where c is dimensionless and at most of order unity. 

Writing our model cut from Eq. G4.2| ) in terms of the variable z, we have 



\ (l-z)(l-z a ) (l-z)(l-z b ) J K ' 

Let f(z) stand for any of the form factors, with corresponding functions <fi(z) from 
Eq. (|2.15|) and P(z) from Eq. (|3.1| ). Consider the difference f(z) = f(z) — g cut (z), and 
let / cut = g cu t<pP. The function fcf>P is thus designed to be analytic on the unit disc. 
We proceed in two steps: First we find a bound for ftp analogous to that in ( |2.14| ); this 



constraint translates into a new bound on the parameters in our expansion. Then we show 
that / cu t is well approximated in the physical region by a polynomial of low degree, so that 
truncating our expansion after only a few terms incurs a very small error. 
By the Minkowski inequality and (|2.14|) we have 



\J d9\U\ 2 j <[J d9\f^j +\J ^|/cut| 2 J < V^a + I^ 2 ), (4.4) 



where z = e ld and 



Y 1-2-n 



I cut = — d6\f cut \ 2 . (4.5) 
271" Jo 



As before, the functions P(z) are unimodular on the unit circle, and so leave the integrals 
unchanged. 7 C ut can be computed using the explicit form for the cut in Eq. (fO[). As a 
specific example, take for cf> in Eq. ( |2.15| ) the case of p — 3 and s = 1, corresponding to 
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the form factor g. The numbers to follow are specific to the case m = mo*, although the 
qualitative results are the same for m = mo- For z a = —0.32 and Z(, = —1.0, we find 

I cut w 2.2 x 10" 3 c 2 . (4.6) 

Thus the bound on / is relaxed only by ~ 5% times c relative to that on /. A realistic 
choice of coefficient c significantly improves this bound, as does a branch point closer to 
the B-D or B-D* threshold z = —1. For example, for z a = —0.5 replace 2.2 x 10 -3 by 
4.0 x 10" 4 in Eq. flO|) . 

Being analytic, f<pP has a Taylor expansion for \z\ < 1. Hence one may write 

^=^(£-"H- (4 - 7) 

Here the coefficients a n are bounded, J2^=o \ a n\ 2 < (1 + lilt) 2 - Moreover, f cut is analytic 
over the physical region. Define the remainder R^ ut through 

oo N 

f cut (z) = = E + R ™t( z ) ■ (4-8) 

n=0 n=0 

The remainder R^ ut is the additional error introduced in the parametrization of f(z) as an 
A^-th order polynomial in z with coefficients a n + b n . These coefficients obey a constraint 
almost identical to Eq. (pOI), because the b n /c are uniformly small: Yln=o l a ™ + ^l 2 — 
(l + llS) 2 (l + En=o Kl) 2 = (1-07) 2 , (1.07) 2 , and (1.08) 2 , for c = 1 and N = 2, 3, and 4, 
respectively. The maximum of R^ ut (z) over the physical region is — 1 x 10 -6 , 2 x 10 -7 , and 
4 x 10 -9 times c for N = 2, 3, and 4, respectively. These figures should be compared with 
the bound on the truncation in the analytic part, J2^=n+i a nZ n < (1 + / ( ! u / t 2 )(0.056) Ar+1 ~ 
2 x 10 -4 , 1 x 10 -5 , and 6 x 10 -7 . Since we expect realistic branch cuts to have c C 1, we 
see that their effect is negligible. 



5. Uncertainties 

Two important statements follow immediately from Eqs. ( |3.2|) and (|3.3|) : 
a ) Each of the various B — > D or B —>■ D* form factors can be accurately described by 
three parameters, one of which is the normalization at zero recoil, with a truncation 
error of order 1%; 
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b ) The fitting parameters obey E™~Q|a n | 2 < B 2 , with B = 1. 

A number of approximations have been made in deriving these results. How do cor- 
rections to these approximations alter the above statements? We answer this question by 
noting that nearly all the corrections we expect to be non-negligible can be taken into 
account by altering the bound to 5 ^ 1. 

An estimate of how much B might change may be made by considering uncertainties 
arising from the following sources: 

1 ) The b and c quark masses, which enter into the one- loop perturbative functions Xji are 

not well established; we take m c /mb = 0.33. This leads to a roughly 5% uncertainty in 
the normalization of <p which, by redefining the a n , is equivalent to a 5% uncertainty 
in the value B bounding the a n . 

2 ) The functions x j a ls° receive perturbative two-loop corrections. Since 4>, and therefore 

the bound, depends only on (Xj) 1//2 , such corrections should lead to no more than a 
15% change in B. 

3 ) The masses of the B* poles were computed from a potential model. Computations 

by various groups typically agree to a fraction of a percent [PH 131- The results from 
two different groups |ll|,^2j give Blaschke factors that agree to 2%. However, P(z) is 
sensitive to the mass of the 3S vector state, which is close to the B-D* threshold, 
and is presented only by fpT2[| . Changing it by 1% results in a 20% change in -P(O); 
P(z)/P(0), however, varies by less than 1%. 

4 ) We argued in the previous section that contributions from multi-particle cuts should 

alter the bound B by less than 8%. 

5 ) In extracting values of | V c b\T{l), we require bounds not on a\ and ai alone, but on 

ai/JF(l) and a^jTiX), and these depend on the zero- recoil normalization This 
normalization is predicted to no worse than 20% accuracy by heavy quark symme- 
try!- 

The uncertainties (1) to (5) are uncorrelated, so to estimate their total effect, we add 
them in quadrature. This leads to a relaxation on our bounds from B = 1 to less than 
B = 1.4. 

This relaxation of the bound increases the truncation error on, for example, fo from 
0.012 to 0.017, still negligibly small given the current experimental accuracy. Even if we 
added the uncertainties linearly, the truncation error would only rise to 0.020. We see 
that statement (a) is extraordinarily robust; it is nearly independent of the size of the 
uncertainties listed above. 



11 



On the other hand, allowing a value of B larger than 1 in statement (b) could in 
principle affect the extraction of \V c b\ significantly due to a larger allowed range for a± and 
a,2- For this reason, it would be useful to pin down the a s corrections and the mass of 
the 35 vector B* state more precisely. However, for the extraction we perform in Sec. 7, 
relaxing the bound from B = 1 to 1.4 turns out to change the results very slightly: The 
central values of | V^t» | and the slope change by no more than a tenth of a standard 
deviation. 

One should also bear in mind that our bounds can be significantly improved by the 
inclusion of more terms than just B-D* or B-D pairs to saturate the bound in Eq. ( |2.8| ). 
Such contributions arise through higher resonances of the current J; if estimated numeri- 
cally, they have an effect equivalent to reducing B. 



6. Heavy Quark Symmetry 

The parametrizations Eq. ( |3.2j ) make no use of heavy quark symmetry. Thus, 1 / m c 
corrections to the extraction of |V C &| from B — > Dlv decays enter only through the nor- 
malization of the form factor f+(z = 0) at zero recoil. This normalization is determined 
by heavy quark symmetry to 0{l/m c ). 

If the individual B — > D*lv form factors /o, Fi, and g are experimentally deter- 
mined in the near future, separate extractions of \V c b\ can be made for each form factor. 
These extractions will depend on heavy quark symmetry only through the normalization 
of form factors at zero recoil. This is useful because the normalization of fo is predicted 



to 0(l/m^)p 



At present, to extract \V c f,\ from the B — > D*W differential width (2^3) in terms of our 
three-parameter descriptions, one must relate fo, a+, and g using heavy quark symmetry. 
In the infinite mass limit, all form factors for B — > D and B — > D* (as well as B — > B) 
are directly proportional to the universal Isgur-Wise function[|J. Consequently, the ratio 
of any two form factors assumes a simple form: 

a+/g = ~, f /a + = -2M 2 r(u + l), and f /g = M 2 r(u + 1), (6.1) 

where u = v -v' is the product of the B with D or D* meson velocities. These ratios admit 
two types of correction, namely power corrections in l/m c , and running and matching 
corrections relating QCD to the heavy quark effective theory. We discuss each of these in 
turn. 
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Because many heavy quark symmetry- violating contributions cancel in the above ra- 
tios, one might expect 1 / m c corrections to be smaller than in, say, the relation between the 
B — > B elastic and B — > D or B — > D* semileptonic form factors. For example, a+/g = — \ 
may be derived using only charm quark spin symmetry, without recourse to bottom-charm 
flavor symmetry; spin symmetry is expected to hold more precisely than the full flavor-spin 
symmetry fill . In addition, the ratio fo/g = M 2 r(uj + 1)[1 + {u — 2)A/2m c ] involves no 
unknown w-dependent functions Hl5"H at (9(l/m c ), but only the constant A = m£,—m c . The 
third ratio is given by the quotient of these two. Choosing two different pairs of the above 
ratios gives two different parametrizations of the decay form factor J-"(uj) conventionally 
defined by 

£ = ^r^ M " mfV^llMu + l) 1 "^/ + (u + 1) 2 ]^V). (6-2) 

In the heavy quark limit, the form factor ^(uj) is simply the Isgur-Wise function times 
QCD corrections (discussed below), and we readily see that Eq. ( |2.3|) reduces to Eq. ( |6.2|) . 
In terms of the parametrizations Eq. ( |3.2|) of g and /o, respectively, 



= d 7 U i \ [^2(0)^2(0)^(2 = 0) + M^F( ai z + a 2 z 2 )} (6.3a) 

P 2 (z)(p2{z) 

[Po(0)M0)Hz = 0) + -±-=(b lZ + b 2 z 2 )}. (6.36) 



The form factor may be expressed as a function T{uS) of velocity transfer by rewriting z 

as 

^+i _ 1 

= • (6-4) 

At zero recoil, J-{ui = 1) = 1 times corrections whose estimates range from 0.89 to 0.99||. 
Relative to this normalization, the parametrization in Eq. ( |6.3a| ) has a 0.5% truncation 
error, while that in Eq. ( |6.3fc| ) has a 1.0% truncation error; see Table 2. 

To the degree that l/m c corrections are negligible, the extracted values of | V c b|jF(o; =1) 
and the slope T'{u = 1) cannot depend on which of the parametrizations (6.3) we use. 
Since l/m c corrections enter differently into each of these parametrizations, the degree 
to which this is true gauges the sensitivity of the extraction to heavy quark symmetry 
violations. 

For a thorough accounting of relations between form factors when using heavy quark 
symmetry, one must also include effects due to the running of the QCD coupling a s and 
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matching between the full theory of QCD and the heavy quark effective theory. The form 
factors are then no longer just trivial factors times the Isgur-Wise function, but now include 
a functional dependence on w, as well as m c , rrib, and the value of a s at these scales. For 
conciseness and definiteness, we adopt the notation of NeubertjTT] to parametrize such 
corrections. The relation between the Isgur-Wise function and the relevant form 

factors then reads 

a+ = -^=(C! + C!r + C!)^ /o = MVr(w + 1)<7& g = jJ^Ctf. (6.5) 

The functions C\, C\ become unity when the strong coupling is switched off, whereas the 
other C"s vanish. In this limit we recover the ratios in Eq. ( |6.1|) . 

Apart from changing the overall normalization of form factors at zero recoil by a 
few percent, the functional dependences in Eq. (|6.5| ) turn out to be rather weak over the 
allowed range for B — > D or B — > D* semileptonic decay {u = 1.0 to 1.5). In particular, 
C\ decreases from 1.136 to 1.011 over this range, but —2a+/g = 0.864 — ► 0.882, and 
fo/gM 2 r(uj + 1) = 0.868 — > 0.884 are nearly constant. In addition, corrections due to 
running between the bottom and charm mass scales cancel out of such ratios. 

Because the undetermined l/m c corrections are just as significant, there is little to 
be gained in incorporating the calculated matching corrections explicitly in our analysis; 
rather, our sensitivity to both l/m c and matching corrections is gauged by comparing 
the extractions of | V^l^l) and the slope •F'(l) by the two parametrizations of Eq. (6.3). 
Compared to the g parametrization, the fo parametrization changes the central values 
of CLEO's | V cb and both CLEO's and ARGUS's ^'(1) by less than a fourth of a 

standard deviation; ARGUS's and ALEPH's | V^l^l), as well as ALEPH's change 
by less than a tenth of a standard deviation {i.e., < 2% for all |V c b|.F(l)). 



7. Results 

Since both parametrizations (6.3) give essentially the same results, we choose the g 
parametrization, which has a smaller truncation error. From the point of view of heavy 
quark symmetry, one should use the fo parametrization, since fo(u = 1) is predicted 
to higher accuracy. Here we are more concerned with exploring the implications of our 
parametrizations. The central values and 68% confidence intervals should be taken as 
indicative; proper inclusion of efficiencies, resolutions, and correlated errors can only be 
done by the experimental groups themselves. 
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Fitting ai/jF(l), and a^jTiX) to experiment yields the results in Table 3. 

For each experiment, we have listed the best fit values for | V^,|.F(1), ai/JF(l), and a^jTiX), 
as well as the resulting slope = 1). The 68% confidence intervals due to statistics 

are included as well. The parametrization ( |6.3a| ) includes the constraint E^°=o l a «J 2 — 1> 
for comparison, we also present the best fit values resulting from an unconstrained fit with 
freely varying a n . 



B 


\V cb \F(l) • 10 3 


ai/^(l) 


a 2 /^(l) 




Expt. 


1 

00 


33.3+JJ 


0.046l° ;° 1 I 
0.181±g;|? 


-1 oo +2 -° 

-3.20±|'| 


-0.89t°; 3 


CLEO 
CLEO 


1 

00 


45.8if£ 9 
49.5^ 


-0.200^;^ 
-0.297^ 


0.981% 
2.59 


9 q+1.2 

-2.8l?; 9 8 


ARGUS 
ARGUS 


1 

00 


3i-5±J:i 
3i.s±J:i 


0.073+° f 3 


1.001% 
1-33±?;1 


-0.65±J;S 
— 0.74±J;| 


ALEPH 
ALEPH 



Table 3. Fit values for |V^,|.F(1), a\/F(l), a2/J-(l), and the zero recoil slope of 
F(oo) from the various experiments, constrained to obey E"^Q|a n | 2 < B. 

The fits allowed by QCD are those with (in particular) \a%\ < 1. The extracted 
values of \V c b\ are in good agreement with a previous extraction^, after accounting for 
differences in definitions and experimental data. We have renormalized the ARGUS data 
to bring their assumed B lifetime and D° — > K~tt + branching ratio into agreement with 



more recent experiments; we use tb = 1.61 psec and ||18|| B(D° — > K it 



4.01%. 



The central values for |U c b|^ r (l) agree surprisingly closely with those of the experimen- 
tal groups themselves. This did not need to be the case, as one can see from the behavior 
of the unconstrained fit. 

The connection between our parameters a±, a<2 and the commonly used expansion in 
(u — 1) is 



^(i) 



1 + 



5.54 



- 1.15 



-7.73-^r + 0.69- ° 2 



+ 1.11 



1' 



+ 



8.19 



1.14 



_02_ 



0.99 



(u,-l) 3 + ... 



(7.1) 

While such an expansion describes the form factor well close to zero recoil, it converges 
poorly over the rest of the kinematic range. Substituting the allowed range of parameters 



E 



n=2 
n=0 



a n \ 2 < 1 gives a truncation error for a quadratic fit in (u>— 1) of 120%; the truncation 
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error of a linear fit is 220%. To be assured of fitting a QCD-allowed form factor at percent- 
level accuracy, a parametrization obeying the same constraints as Eq. (6.3) must be used. 

Plotted in Fig. 1 are the constrained and unconstrained fits to the CLEO0 data. 
Both fits match the data well; the chi-squares per degree of freedom are x 2 / dof = 0.65 
and 0.50, respectively. The CLEO group extracts • 10 3 = 35.1 ± 1.9 (stat) and 

a slope ^'(1) = —0.84 ± 0.13 using a linear fit, in close agreement with our bounded fit. 
The unbounded fit serves as an illustration of a parametrization which gives a markedly 
different best fit; the central value of |V c f,|jF(l) differs by 5% from the linear result, while 
the slope is in violation of the Bjorken boundpj, .F'(l) < -1/4. By Eq. flPD , this 
uncontrained fit is ruled out by QCD. 



45 



40 



CLEO 




Figure 1. Fit of CLEO data to our parametrization, Eq. ( |6.3a ). The solid line shows the 
result of imposing the QCD-derived constraint ^ 



n=2 
n=0 



< 1 on the parametrization. 



The dot-dash line shows the corresponding unconstrained fit. 



For ALEPHQ the constrained and unconstrained fits overlay each other quite closely 
(Fig. 2). A linear fit by the ALEPH group gives | V c6 |.F(l) ■ 10 3 = 31.4 ± 2.3 (stat) and a 
slope ^'(1) = —0.39 ±0.21, in good agreement with the results of our constrained fit. The 
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confidence intervals in Table 3 are somewhat larger for ALEPH than might be expected 
because of the smallness of the minimum % 2 : Both the bounded and unbounded fits have 
X 2 / dof = 0.37, so a larger range of fit parameters fall within the 68% confidence limits in 
either case. 

45 
40 
35 

>^ 

30 
25 
20 

1 1.1 1.2 1.3 1.4 1.5 

CO 

Figure 2. Fit of ALEPH data to our parametrization; see Figure 1 caption for details. 

The constrained and unconstrained fits to ARGUS data differ mainly near zero 
recoil, with comparable values x 2 1 dof = 0.70 and 0.67, respectively (Fig. 3). The ARGUS 
group used several parametrizations, which yielded central values of |U c b|jF(l) ■ 10 3 from 
39 to 46. Their linear fit gave |U cb |J^(l) ■ 10 3 = 39 ± 4 and = -1.17 ± 0.11, in some 

contrast to our central values. 
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Figure 3. Fit of ARGUS data to our parametrization; see Figure 1 caption for details. 

Although the large statistical uncertainty in a 2 /^(l) precludes its determination at 
present, we can make a definite prediction for the future: The central value of aijTiV) 
must increase from CLEO's present (unconstrained fit) number to fall inside our bounds. 
Taking the theoretical estimates of Sec. 5 into account, we predict \ailT{\) \ < 1-4- 



8. Conclusions 

Dispersion relation techniques and the use of analyticity properties of hadronic form 
factors as functions of their kinematic variables provide a valuable window into the realm of 
nonperturbative physics. Using these methods, one can obtain useful bounds on quantities 
of interest, in this case the form factors in the semileptonic decays B — > DW and B — > D*W. 

These bounds may be transformed into parametrizations of the four experimentally 
accessible form factors relevant to B — > DW and B — > D*W. Given the continuing experi- 
mental scrutiny devoted to these decays, these form factors will likely be measured in the 
foreseeable future. 
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Our derivation of these parametrizations relied on dispersion relations, crossing sym- 
metry, and a perturbative QCD calculation performed at a scale mg+m]). . The derivation 
improves on an earlier work|7j in that no use of heavy quark symmetry was made. The 
various uncertainties involved in the derivation, such as perturbative corrections and un- 
certainties in quark masses, were estimated, and shown to be unimportant. This includes 
effects from branch cuts in the form factors due to non-resonant contributions. 

The result is a three-parameter description of each of the form factors /o, g, and 
/+ accurate over the entire physical kinematic range to better than 2%. The value of 
one of the parameters, the normalization at zero recoil is predicted by heavy quark 

symmetry. The other two parameters a\,ai are constrained by |ai| 2 + 1 0,2 1 2 < B 2 1 with a 
leading-order result B = 1. A very conservative estimate of corrections to our results leads 
us to conclude that to all orders, the bound obeys B < 1.4. The 2% or better accuracy of 
the three-parameter description applies for any B < 1.4. The three parameter fit to the 
form factor F\ is less accurate; for B < 1.4 we find a bound of 8% on its relative error. 

We emphasize that we have determined strict upper bounds on the truncation errors. 
The truncation errors may be significantly smaller. The strict inequality ( |2.8| ) can be 
improved by including the contributions of other intermediate states; our use of Blaschke 
factors, Eq. ( |3.1|) , amounts to assuming the largest possible uncertainty from the residues of 
poles in the form factors; and the bounds on the parameters for each form factor, Eq. ( |3.3|) , 
are actually correlated, as in Eqs. (|3.4|) and (|3.5[) . 

As an application of our results, the individual form factors in B — > D*W were com- 
bined using heavy quark symmetries in order to obtain a single parametrization of the 
differential cross-section dF/dq 2 , which was then fit to data. This was necessary because 
the best data currently available sums over D* mesons in all polarization states and thus 
involves more than one form factor. However, to 0(l/m c ), our results depend only on 
charm quark spin symmetry and the constant A = mo — m c , and are therefore expected 
to be more reliable than those using the full flavor-spin heavy quark symmetry. We obtain 
values for a three-parameter (| V^l^l), ai/jF(l), and 02/^(1)) fit to the single form factor 
T{y ■ v') describing B — > D*W that is free of the theoretical errors inherent in choosing 
a parametrization for extrapolating the data to zero recoil. We again emphasize that, 
although heavy quark spin symmetries were used in obtaining values for | V c b\T{l), this 
is a limitation imposed by the currently available data that will be lifted when better 
measurements of the individual form factors become available. 
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The intensive experimental effort focused on semileptonic B — > DW and B — > D*W 
decays will result in increasingly precise measurements of the rate and form factors. Our 
descriptions of these form factors are remarkably insensitive to theoretical uncertainties, 
and are accurate over the physical kinematic range to better than 2%; as such, they should 
be useful ingredients in studying the nonperturbative physics of semileptonic B decays. 
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